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We investigate the graphene-impurity interaction problem by combining experimental - scan- 
ning tunneling microscopy (STM) and spectroscopy (STS) - and theoretical - Anderson impurity 
model and density functional theory (DFT) calculations - techniques. We use graphene on the 
SiC(0001)(2 X 2)c reconstruction as a model system. The SiC substrate reconstruction is based 
on silicon adatoms. Graphene mainly interacts with the dangling bonds of these adatoms which 
act as impurities. Graphene grown on SiC(0001)(2 x 2)c shows domains with various orientations 
relative to the substrate so that very different local graphene/Si adatom stacking configurations can 
be probed on a given grain. The position and width of the adatom (impurity) state can be analyzed 
by STM/STS and related to its local environment owing to the high bias electronic transparency 
of graphene. The experimental results are compared to Anderson's model predictions and comple- 
mented by DFT calculations for some specific local environments. We conclude that the adatom 
resonance shows a smaller width and a larger shift toward the Dirac point for an adatom at the 
center of a graphene hexagon than for an adatom just on top of a C graphene atom. 

PACS numbers: 81.05.ue, 68.35.Dv, 68.37.Ef, 31.15.A- 



I. INTRODUCTION 

Graphene is a one-atom thick crystal composed of car- 
bon atoms arranged on a honeycomb lattice. From this 
crystallographic structure with two equivalent carbon 
atoms - labeled A and B- per unit cell arises a unique 
low-energy electronic structure. It is characterized by 
a linear and isotropic dispersion relation around the K 
(K') corner of the Brillouin zone [Ul^]. This structure 
is referred to as Dirac cones and the low-energy excita- 
tions in graphene behave like a 2D gas of massless Dirac 
fermions [Hill]- Among other properties, high carrier mo- 
bilities that have been measured in this system make it 
promising for applications in electronics 

Except for very specific cases, graphene is in interac- 
tion with its environment. Substrate is known to af- 
fect the graphene properties [21 [Z] as well as intrinsic 
or intentionally created defects [8l412j. Focusing on the 
latter, graphene is expected to show unusual behavior 
when interacting with defects. Going more into details, 
theory predicts an anomalously small broadening of the 
adatom level. When considering magnetic impurities, 
this property makes magnetic moment formation eas- 
ier |13j . Moreover, the adatom adsorption site may in- 
fluence the physics of the Kondo effect [TH [TS]. The 
adsorption-site dependency of the impurity level broad- 
ening and shifting has been computed |16H19j . providing 
an impurity adsorption-site signature that can be mea- 
sured by STM. Yet no experimental observation of this 
effect has been given so far. 

In this study, graphene is obtained by graphitization 
of SiC [201 in]- We consider the graphene-SiC(OOOT) 
interaction and we show how our results are relevant 
to the investigation of the graphene-impurity problem. 
When graphene is grown on SiC(OOOl) under ultra high 
vacuum, two SiC surface reconstructions (3 x 3) and 



(2 X 2)c are found at the interface [52]. In case of 
the (3 X 3) interface, the graphene-substrate interac- 
tion is very weak [52 [33] while it is stronger for the 
(2 X 2)c [2211211125]. The n states are however always 
detected around the Fermi level in this system ^22; . Here 
we focus on the SiC(0001)(2 x 2)c interface. The atomic 
structure of the (2 x 2)c reconstruction consists in one 
Si adatom and one C restatom per unit cell, displaying 
an empty/filled DB respectively. Our study deals with 
the hybridization between the Si adatom level - which 
plays the role of a defect - and the graphene states. 
Graphene islands being rotated with respect to the sub- 
strate, spatial variations of the graphene/adatom stack- 
ing are observed. This allows us to experimentally in- 
vestigate the stacking dependence of the interaction for 
a given impurity level, which is hardly achievable for de- 
posited adatoms since they usually have a well defined 
adsoption site [23 123 ■ The specificity of STM imaging 
of graphene on SiC, with high-bias measurements giving 
access to the interface states and low-bias measurements 
to the graphene states [2S], allows us to investigate the 
impact of the interaction on both the adatom and the 
graphene states. 

In section |ll| we provide information on the sample 
fabrication, on STM measurements and on the ab ini- 
tio calculation parameters. Section III begins with an 
introductive paragraph, providing more details on how 
graphene on SiC(000T)(2 x 2)c allows us to study ex- 
perimentally the graphene-impurity interaction problem. 
We then present our experimental and theoretical results. 
Low bias STM images display stacking dependent per- 
turbations of the graphene local density of states at low 
energy, the stronger (weaker) perturbations appearing 
for top (hollow) graphene/adatom stacking. High bias 
constant current STM images together with tunneling 
spectroscopy reveal stacking dependent modifications of 
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the adatom resonance in terms of broadening and shift- 
ing. We then show that the Anderson single-impurity 
model accounts for the general trends we observed ex- 
perimentally. Finally, we complement our study of the 
graphene/adatom interaction with ab initio calculations 
for the whole graphene/reconstruction system. Keeping 
in mind the increased system complexity and the fact 
that we consider one specific geometry, we concentrate on 
the general trends of the computed electronic structure. 
We identify a stacking-dependent impact of graphene on 
the adatom resonance similar to the one derived from 
the simple Anderson model. Regarding graphene states, 
maps of integrated DOS show features similar to the ones 
observed on low-bias STM images. 

II. EXPERIMENT AND CALCULATION 
DETAILS 

A. Experimental aspects 

The sample preparation and characterization were per- 
formed under ultrahigh vacuum. We followed the SiC 
graphitization procedure presented previously |22| . We 
used n doped 6i?-SiC(0001) samples priorly cleaned by 
annealing under a Si flux at 850° C. Annealing steps at 
increasing temperature up to 1100° C allow us to get a 
graphene coverage of less than a monolayer. Typical 
LEED patterns show SiC(3 x 3) and SiC(2 x 2) spots 
and a ring-shaped graphitic signal with modulated in- 
tensity [HIMllMllSn]- 

The STM and STS measurements were conducted 
at room temperature with mechanically cut Ptir tips. 
On the surface of the samples, we find domains show- 
ing the bare SiC(0001)(3 x 3) reconstruction, graphene 
monolayer islands on the SiC(0001)(3 x 3) reconstruction 
(G/S X 3) and on the SiC(000T)(2 x 2)c reconstruction 
{G/2 X 2) and also few multilayer islands, similarly to 
previous studies [22l [23l [30] . In this paper, we will fo- 
cus on the G/2 X 2 islands using constant current im- 
ages and tunneling spectroscopy in the current-imaging- 
tunneling spectroscopy (CITS) mode. In this latter case, 
a topographic image is acquired at the stabilization con- 
ditions {Vstab,It)- At each point, the current feedback 
loop is opened while the sample-surface distance (defined 
by [Vstabi It)) is maintained constant, in order to record 
an I(V) curve. The dl/dV curve is then numerically cal- 
culated. 



B. ab initio calculation framework 

Calculations are carried out using the VASP code [3T] , 
which is based on the density- functional theory (DFT). 
We use the generalized gradient approximation |32j to- 
gether with a plane wave basis and ultrasoft pseudopo- 
tentials [M^. The 4H-SiC substrate is modeled by a slab 
containing 4 SiC bilayers. The backside of the slab is 



passivated by hydrogen. An empty space of 8 A sepa- 
rates the graphene layer from the next SiC slab. We use 
a plane wave basis cutoff of 211 eV. The ultrasoft pseu- 
dopotentials have been extensively tested [211[311[3S]. In- 
tegration over the Brillouin zone is carried out within the 
Monckhorst-Pack scheme, using a6x6x5tol8xl8xl 
grid. All the structures were fully converged, with resid- 
ual forces smaller than 0.015 eV/A. The local density 
of states (LDOS) is calculated by projecting the Kohn- 
Sham wave functions onto site-centered spherical har- 
monic functions within a sphere of radius r = 1.65 A 
centered on the considered site. 



III. RESULTS AND DISCUSSION 

A. From graphene on 6H-SiC(000T)(2 x 2)c to the 
impurity problem 

The system we study experimentally is composed of a 
monolayer of graphene on SiC(0001)(2 x 2)c- The atomic 
structure of the bare substrate reconstruction is now well 
established and is given in figure [TJa). As it was first pro- 
posed from quantitative LEED measurements and STM 
by Seubert et al.[3B], three out of the 4 dangling bonds 
of the (2 X 2) bulk-truncated cell are saturated by a Si 
adatom, a priori leaving two half-filled dangling bonds 
per cell, one on the Si adatom and one on the unsat- 
urated C atom. Later, ab initio calculations [24] along 
with STM [25l [35] revealed a charge transfer from the Si 
adatom to the unsaturated fourth C atom, leaving one 
empty dangling bond on the Si adatom and one filled 
dangling bond on the C atom, which we call a restatom. 

In this work, we address the question of the hybridiza- 
tion between graphene tt states and states from the sub- 
strate reconstruction. A schematic side view of the 
system is given in figure [l] (b). We mainly focus on 
the graphene/adatom interaction since previous ab ini- 
tio calculations of the band structure of graphene on 
SiC(0001)(2 X 2)c have shown that the interaction be- 
tween graphene and the C restatoms is much weaker than 
with the Si adatom - calculations were done either us- 
ing a GGA ^24J or a GGA(D) [37] exchange-correlation 
functional, the latter one taking van der Waals interac- 
tions empirically into account. This statement is rather 
intuitive since as depicted in figure [l](b). Si adatoms are 
much closer to the graphene layer than C restatoms, cal- 
culations giving an adatom/restatom height difference of 
0.9 A [M". 

Another characteristic of the system which is crucial 
for our study is that on the C face of SiC, graphene 
grows with a strong rotational disorder under UHV and 
typical LEED patterns show the graphene signal as a 
ring with modulated intensity (see EPAPS figure SI 
(a)). Figure [l] (c) gives a schematic representation of 
the graphene/adatom stacking for a rotation angle of 
a = 16°. The common pseudocell (given in black) con- 
tains various graphene/adatom stacking configurations. 
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The corners of the unit cell display hollow stacking i.e. 
the adatoms sit under the center of a graphene hexagon. 
Within the cell, we identify top (shifted top) stacking 
i.e. the adatoms sit under (close to) a graphene atom 
or bridge stacking i.e. the adatoms sit under a C-C 
bond. Unlike in a graphene-adsorbate system, the po- 
sition of the reconstruction adatom with respect to the 
graphene lattice is imposed by the interface geometry 
so that we can investigate the graphene-adatom interac- 
tion using room temperature STM measurements. More- 
over, adsorbates usually sit on well defined (top, bridge or 
hollow) energetically favored surface sites [IHl HZj which 
would make the investigation of the stacking-dependence 
of the interaction hardly possible. In case of (7/2 x 2, 
spatial variations of the graphene/adatom stacking are 
imposed by the system and easily identifiable thanks to 
the transparency of graphene in high bias STM images, 
as will be discussed in section IIIIBI in connection with 
figure [2j Note finally that the adatom states can be con- 
sidered as independent in a first approximation because 
they are separated by 6.14 A and the graphene/adatom 
interaction is still sufficiently weak to prevent a coupling 
through graphene states. Hence, the system provides the 
opportunity to investigate how the interaction between 
graphene and a localized impurity state depends on the 
local configuration without being limited to stable adsor- 
bate positions. 



B. Scanning tunneling microscopy and 
spectroscopy results 

Figure [2] (a) and (b) present dual-bias STM images 
of a G/2 X 2 island {Vs = ±1.5 V). At such high voh- 
age bias, graphene appears transparent and the images 
show the reconstruction states, namely with one protru- 
sion per unit cell. At negative (positive) sample bias Vs, 
we probe filled (empty) states. The reference cell (dia- 
mond cell) reveals the spatial shift between filled/empty 
states which is consistent with charge transfer between 
the reconstruction adatom and restatom predicted by ab 
initio calculations of the electronic structure of the bare 
SiC(000T)(2 X 2)c reconstruction |24|- Empty (filled) 
states are ascribed to states located on the Si adatom 
(C restatom) from the mentioned ab initio calculations. 
The electronic structure of the substrate reconstruction 
is only weakly affected by the graphene layer, indicating 
a small or moderate graphene-adatom interaction. 

Figure [2] (c) and (d) are dual-bias STM images of a 
G/2 X 2 island at high Vs > and low Vs, respec- 
tively showing the reconstruction adatom states and the 
graphene atomically resolved LDOS near the Fermi level 
{Ep). On both images, we identify a superstructure of 
periodicity P = 3.0 ± 0.3 nm. All the G/2 x 2 islands 
observed show such superstructures, however with vari- 
ous periodicities in the nm range. We interpret them as 
moire patterns that arise from the composition of the 
graphene and the reconstruction lattices stacked with 




FIG. 1: (Color online) (a) Schematic representation of 
the SiC(0001)(2 x 2)c surface reconstruction. The unit cell 
(dashed diamond cell) contains one Si adatom with an empty 
DB (dashed line) and one C restatom with a filled DB (dot- 
ted line), (b) Schematic side view of the system G/2 x 2 with 
the Pz graphene orbitals and the empty (filled) DB on the 
reconstruction adatom (restatom). (c) Schematic top view 
of the graphene/adatom stacking in case of a rotation angle 
a = 16°. The common pseudocell (black solid line) contains 
various stacking configurations. 



rotation. The geometrical proof of this interpretation, 
based on the measured relation between the moire pe- 
riodicity and the graphene/substrate rotation angle, is 
given in the supplementary information and is analogous 
to the one we applied in [23| to G/3 x 3 islands. Thus, 
the moire pattern (periodicity P — 3.0 ± 0.3 nm, stack- 
ing angle a = 16 ± 1°) is visible on both the adatom 
and the graphene lattice (pseudocell in solid line). On 
the adatom sites, it appears as a periodic spatial modu- 
lation of their apparent height and on the graphene lat- 
tice as a periodic repetition of perturbed/unperturbed 
LDOS zones. The graphene perturbations consist in a 
LDOS deficit on graphene atoms ("switched off" atoms) 
distributed with a periodicity similar to the one of the 
underlying reconstruction [22]. We recall that, from a ge- 
ometrical point of view, the moire pattern is associated 
to the quasi periodicity of the graphene/reconstruction 
stacking. Hence, the perturbations observed on both 
the adatom and graphene lattice being correlated to the 
moire pattern, the graphene substrate/interaction should 
vary with the local stacking. 

To further study the stacking dependence of 
the graphene/adatom interaction, we can relate the 
graphene/adatom stacking and the perturbations using 
the fact that figure [2] (c) and[2](d) are dual bias images. 
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We can thus identify directly the stacking which corre- 
sponds to the "switched off" graphene atoms. Figure 
[l] (c) in the previous section corresponds to an iUustra- 
tion of the graphene/adatom stacking obtained this way, 
using figure [2] (c-d) . On the corners of the pseudocell 
(in sohd fine), the adatoms sit under graphene hoUow 
sites. This particular stacking corresponds to small re- 
gions (sa 1 nm wide) showing honeycomb contrast typical 
of graphene on the low bias STM image. Within the two 
halves of the pseudo cell, adatoms sit under one of the 
two graphene sites in a top, shifted top or bridge con- 
figuration. Those configurations lead to reduced LDOS 
("switched off" atoms) on the concerned graphene sites, 
in agreement with our previous ab initio calculations [21] . 
Thus, from the point of view of graphene states, the im- 
pact of the graphene/adatom interaction is very weak 
in the hollow stacking configuration while it is stronger 
in the top, shifted top and bridge configurations. In the 
following, we will concentrate on the effect of this interac- 
tion on the adatom states and analyze the moire pattern 
on high bias images. 




b) 



d) 



FIG. 2: (Color online) (a) and (b) 3 x 3 nm^ Dual-bias STM 
images of a G/2 x 2 island. At high bias, images show the 
substrate reconstruction states (solid line unit cell), (a) ((b)) 
exhibits the filled (empty) states associated to the restatom 
(adatom). (c) and (d) 6x6 nm'^ Dual-bias STM images of 
the same Gjl x 2 island as in (a-b). (c) exhibits the adatom 
empty states while (d) exhibits the graphene overlayer LDOS 
near Ep. A moire pattern (solid line cell) is visible on both 
images (periodicity P = 3.0±0.3 nm, measured stacking angle 
a = 16± 1°). 

STM images contain information on the topography 
and the electronic density of states of the sample surface. 
In order to decorrelate those two contributions and inves- 
tigate the nature of the moire pattern that appears when 
imaging the adatom states, we first focus on the depen- 
dence of constant current STM images on the tunneling 
bias. Figure [3] presents a series of STM images taken 
at various sample biases on a G/2 x 2 island. System- 



atic dual-bias image acquisition at the reference tunnel- 
ing bias Vs = -1-1.5 V (see figure |3] (a)) ensures that im- 
ages are taken at the same spot of the island and that no 
tip changes occurred between images. The island shows 
an incommensurate moire pattern (black pseudocell) of 
pseudoperiodicity P — 3.3 ±0.3 nm. We measure the ro- 
tation angle a = 17 ± 1 ° between the graphene and the 
(2 X 2) lattices using the low bias image (Vs = 10 mV) 
in figure |3] (b) . 

The low bias STM image in figure [3](b) {Vs = 10 mV) 
furthermore allows us to identify the graphene/adatom 
stacking and to locate the typical zones of stronger and 
weaker perturbations of the graphene electronic structure 
that we presented in connection with figure [2](d). The 
moire pseudocell is represented in black. On the whole se- 
ries of images the corners of the cell correspond to hollow 
stacking and to weaker graphene perturbation. Within 
the cell, we identify top, shifted top and bridge stacking 
which correspond to stronger graphene perturbation. 

The series of high positive bias images showing the 
adatom states (figures [s] (c-e)) exhibits two inversions 
in the moire pattern contrast between V5 = -1-0.5 V 
and Vs — +2.0 V. This behavior indicates that the 
moire pattern is related to a spatial modulation of the 
adatom LDOS rather than to a topographic modula- 
tion. This aspect will be further investigated through 
scanning tunneling spectroscopy. As a side remark, one 
can note the asymmetry between the two halves of the 
moire pseudocell at Vs = -1-1.0 V, the top half appearing 
darker than the lower one. This differentiation is unlikely 
to arise from the type of graphene sublattice impacted 
since they are equivalent in LDOS in ideal graphene. A 
schematic representation of the graphene/reconstruction 
stacking given in figure |3] (f) helps to find an explana- 
tion to this observation. Both halves of the pseudocell 
show top graphene/adatom stacking at their center, the 
difference between them lies in the graphene/restatom 
stacking. The upper (lower) half shows a top (hollow) 
graphene/restatom stacking. A weak graphene/restatom 
interaction could thus explain the differentiation between 
the two halves of the moire pseudocell. This assump- 
tion is further supported by the differentiation of the re- 
statoms in correlation to the moire pattern revealed by 
the STM image at V^s = -1-0 V (filled states) in figure 
[3] (g) . Finally, this interpretation is consistent with the 
residual graphene/restatom interaction obtained in re- 
cent GGA(D) ab initio calculations for this system j37j . 
which authors show as responsible for n-type doping of 
the graphene layer. 

In the following, we focus on the much stronger 
graphene/adatom interaction and complement data from 
figure |3] with scanning tunneling spectroscopy measure- 
ments. The constant current image contrast is generally 
interpreted as a map of isovalues of the sample LDOS 
ps{{x,y),Ep) for low Vs or of ps{{x,y),E)T{E,Vs.z) 
integrated from Ep to Ep -f eVs for higher tunneling 
biases, with T(E ^Ys-,z) the tunnel barrier transmission 
coefficient. In a spectroscopy measurement, the dl/dVs 



5 




FIG. 3: (Color online) Tunneling bias dependance of the moire pattern (pseudocell in solid line) contrast, 7x7 nm^ constant 
current mode images, (a) STM image at the reference tunneling bias Vs = +1.5 V, (b-e) and (g) being extracted from 
dual-bias images in order to check the structure positions, (b) Low-bias image showing the unperturbed (corner of the moire 
pseudocell) and pertubed (within the moire pseudocell) graphene zones, (c-e) Evolution of the moire contrast with Vs on 
adatom empty states, (d) shows an asymmetry contrast between the top and bottom halves of the moire pseudocell. (f) 
Schematic representation of the graphene/reconstruction stacking deduced from (a), (b) and (g). Honeycomb contrast on (b) 
corresponds to hollow graphene/adatom and top graphene/restatom stacking. Perturbed graphene regions correspond to top 
graphene/adatom stacking with top (hollow) graphene/restatom stacking in the top (bottom) half of the moire pseudocell. (g) 
Moire contrast on restatom states (filled states). 



is considered as a probe of ps{{x, y), Ep + Vs), for slowly Figure |4] shows data extracted from a CITS measure- 
varying T{E,Vst z). ment on the same G/2 x 2 island as in figure [s] The 
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corresponding 7x7 nm^ topographic image at stabiliza- 
tion parameters {Vstab — +2.2 V, It = 0.6 nA) is given in 
figure|4](a). The dl/dV map at Vs = +1.0 V in figure|4] 
(b) clearly reveals a modulation of the surface LDOS in 
correlation with the moire pattern. Figure |4] (c) and (d) 
show I{V) and dl/dV curves acquired within the three 
different stacking domains (indicated by crosses on figure 
|4](a-b)) identified using data in figure ^ namely hollow 
graphene/adatom (red), top graphene/adatom and hol- 
low graphene/restatom (blue) and top graphene/adatom 
and top graphene/restatom (black). In each spectrum, a 
clear feature appears at Vs = +1.35 V, Vs — +1.52 V 
and Vs = +1.61 V respectively [35], which we ascribe 
to the adatom resonance. This establishes that the 
stacking-dependent graphene/adatom interaction shifts 
the adatom resonance energy. The adatom in top config- 
uration shows the higher resonance energy, which corre- 
sponds to configuration leading to the stronger impact on 
graphene low-energy electronic structure (figure [3| (b)). 
Now going back to the data presented in figure [3] the 
shifts in the resonance energy do explain the moire con- 
trast inversion between Vs = +1.0 V and Vs = +2.0 V 
constant current images of figures |3](d) and (e) but does 
not account for the lower-energy contrast inversion be- 
tween figures [3] (c) and (d). This lower-energy contrast 
inversion could however be interpreted as arising from a 
stacking dependent resonance broadening along with the 
stacking dependent shift in energy. The hollow adatom 
resonance would thus be lower in energy but narrower 
than the top adatom resonance, allowing us to probe 
those latter states at energies < 0.5 eV, like in figure 
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The existence of those low-energy states is clear from 
constant current images but hardly evidenced by CITS 
data. An explanation to this comes from the fact that 
each spectrum is acquired at a fixed tip-surface dis- 
tance which is determined by the stabilization parame- 
ters (Vstab, It) while during constant current imaging, the 
tip-sample distance is adjusted in order to keep the tun- 
neling current fixed. As one can infer from the dramatic 
current drop with Vs in figure |4](c), during the acquisi- 
tion of the spectra presented in figure |4] the tip-surface 
distance was relatively high with respect to the one dur- 
ing a constant current image at {Vs = 0.5 V, It = 0.2 nA) 
presented in figure [s] (c). Z{V) curves on this surface ac- 
tually indicate that the tip moves toward the surface by 
« 3 A between Vs = +2.2 V and Vs = +0.5 V (see 
EPAPS). Thus, low bias states are hardly probed in the 
CITS measurements of figure |4] 



In the following section, we will use a simple theo- 
retical model in order to check if a stacking dependent 
graphene/adatom interaction can indeed lead to varia- 
tions in the shifting and broadening of the adatom level. 



FIG. 4: (Color online) CITS data on the same island as 
in figure (sl (a) 7 x 7 nm'^ constant current image at the 
CITS stabilization parameters (moire pseudocell in black), 
(b) Conductance image revealing the moire pattern. (c) 
I(V) curves spatially averaged over 2.5 nm^ (20 spectra), on 
the three different graphene /reconstruction stacking regions 
(see crosses in (a) and (b): graphene/adatom: hollow (red) 
graphene/adatom top and graphene/restatom hollow (blue) 
graphene/restatom top (black), (d) Corresponding conduc- 
tance curves. Each curve shows a peaked structure ascribed 
to the adatom band, at Vs = -1-1.35 V (red), Vs = +1.52 V 
(blue) and Vs = +1.61 V (black). 



C. The Anderson single-impurity model 

In this section, we follow the calculation procedure pre- 
sented in [131 1161 and apply Anderson's description of a 
localized state in the continuum [39 to our system i.e. 
graphene on top of an adatom in the top and hollow 
stacking configuration (see figure [s] (a)). In this frame- 
work, we model the adatom as an isolated localized state 
of energy ojaa and thus neglect the presumably low in- 
teraction between adatoms because of the large adatom- 
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adatom distance -6.14 A- as well as the small impact of 
the restatoms located further away from the graphene 
layer. The spin degree of freedom will also be omitted 
since we are not dealing with magnetic defects (i.e. we 
consider the non interacting Anderson model [40|). Con- 
sequently, the total Hamiltonian of the system is given 
by: 

H=Ug+ Haa + (1) 

Hg is the tight-binding Hamiltonian for isolated graphene 
in the nearest neighbor approximation: 

•Hg = -t^a^(Ri)6(Rj) +^^.0. (2) 

a, b are the fermionic operators for each A and B 
graphene sublattice. Rij are the atomic positions in the 
direct space, t « 2.8 eV is the hopping integral between 
two nearest neighbors. In the momentum space, we get: 

'^9^-tJ2 [^(k)«k^k + <^(k)*6l,ak] (3) 

k 

with 0(k) = Etie'*'^-- Si ^ X, 62 ^ + and 
63 = — are the vectors of the graphene lattice 
connecting a B-type atom to its three nearest neighbors 
(see figure [5] (a)). Haa is the Hamiltonian of the unper- 
turbed adatom orbital of energy ljaa with the fermionic 
operator c : 

Uaa = ^aac'c (4) 

T-lint is the Hamiltonian that describes the coupling be- 
tween the adatom s-like orbital -the adatom s — Pz 
atomic orbital showing axial symmetry- and the graphene 
states continuum. Hint depends on the graphene/adatom 
stacking (see table |l]). For the top stacking CHmt = ^t), 
hopping between the adatom and the graphene sites up 
to the second neighbors is considered. For the hollow 
stacking {Hint — Hh), hopping between the adatom and 
the six nearest graphene neighbors is considered. 





+aV Ek ['^(k)6l. + 0(k)*at ]c + H.c. 


Ht 


+y T,k[bic + «0(k)4c + H.c] 



TABLE I: Graphene/adatom interaction Hamiltonian for the 
hollow {'Hh) and top ("Ht) stacking configurations. 



adatom orbitals with respect to the graphene lattice and 
the local graphene/adatom interaction is thus likely to 
involve more than one graphene atom. Also, this cor- 
responds to the same spatial interaction cut off in the 
hollow and top configuration (see figure [s] (a )). This re- 
finement of the calculation presented in [13|, [T6] is nec- 
essary to treat the graphene adatom interaction on an 
equal footing for the two stacking configurations we con- 
sider here. 

We calculate the localized state retarded Green's func- 
tion Gcc{t) = - (T[c(t)c(0)^]) (with T the time ordering 
operator) whose expression is: 

G^^{u)^[Lu-LUAA~^Mr' (5) 

with Scc('^) the localized state self-energy arising from 
the hybridization to the graphene continuum. The hollow 
(E^(a;)) and top (E^^(cj)) configuration self-energies are 
given in table[ll] as functions of matrix components of the 
bare graphene retarded Green's function G^y'^{uj) with 
x,y = a,b: 

^xvM"^) uj^ -t-^\(l)(k)\-2 + iO+sign{uj) ^ ' 

where cr*' is the identity matrix and {j — 1, 2) the real 
and imaginary off-diagonal Pauli matrices. 





a'V' Ek[<^(k)ea.k(^) + 0(k)*ei,,k(a;)] m 
with 9,,^{cj) = 0(k)GS^k(c^) + </>(k)*G»f,k(a;) 




V' Ek [(1 + c^'m)\')Gt^M(^) + 2a0(k)GL^k('^)] 



TABLE IL Hollow (£^(0;)) and top (Eif^(a;)) configura- 
tion self-energies as functions of the bare graphene retarded 
Green's function coefficients GU^.k('^)- 



We apply the linear-band approximation which is valid 
for sufficiently low energy and express the energy as 
£±(fc) = ivpk with vf « 10^ m/s around the K(K') 
point and chose a cut off energy D k, 7 eV in order to 
match the number of states in the first Brillouin zone. 
We then get the localized state self-energy expressions in 
table Uni 



In these expressions, V is the coupling energy between 
the adatom state and the first-neighbor graphene site in 
the top configuration. We introduce the a < 1 parame- 
ter in order to take into account the distance dependence 
of the coupling energy. aV corresponds to the hopping 
energy between the adatom and the first graphene neigh- 
bors in the hollow configuration or the second neighbors 
in the top configuration (see figure [s] (a)). Taking the 
second neighbors into account in the top configuration 
is motivated by the fact that STM images display large 
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-Lo^f^[D^ +LO^ln\l- ^\]-2tAa^^-^e{D-\Lu\) 


[1311161 


Ej,H, a = 


-^ln\l-^\-iA\co\e{D-\c^\) 






[a^D'' AtaD + 2atLjln + {a^uj^ + t^)ln\l |] 
~i^[a^\Lu\'^ - 2atLU^ + t^\uj\]e(D - \uj\) 





TABLE III: Hollow (E^(a;)) and top (T.'^ci^^)) configuration 
localized state self-energies within the linear-band approxima- 
tion. A = 7r(]^)^ is the dimensionless hybridization parame- 
ter. 



We are interested in the DOS associated to the adatom, 
which in this model reads : Pcc(w) = —^Im[Gf^{uj)]. The 
adatom DOS can be expressed in terms of imaginary and 
real part of the self energy as follows: 

From equation [7j the shift in energy of the adatom 
band due to the coupling to the graphene states is related 
to the real part of the self energy. The renormalized level 
energy ujaa is given by : 

OJAA - UJAA ^ Re[T,{u!AA)] (8) 

Figure [s] (b) represents the calculated interacting 
level energy in both hollow {(jJaa) ^.nd top (ujaa) 
configurations for loaa = 1 eV as a function of a. 
For a > 0.3 , to a A ^aa which corresponds to 
the experimental situation. More generally, by solv- 
ing Re[T,^ {ujaa)] = Re[T,'^ {ujaa)], we get that for 
1 eV < UJAA ^ 2 eV, the switch in the order of the 
resonances occurs for 0.2 < a < 0.3. Thus, our model 
is consistent with the experimental results when a is 
sufficiently large. Actually, considering a ~ 1 is a 
reasonable choice since STM images (figures |2][3| show 
orbitals of axial symmetry on the Si adatoms that 
are large with respect to the graphene lattice - they 
cover approximately one graphene hexagon. Consis- 
tently, ab initio calculated cross sections of (empty 
states) integrated DOS (figure lb in reference [5S]) 
show an orbital of axial symmetry with respect to the 
^^(^[OOOl] axis on the Si adatom, which is large with 
respect to the graphene cell. An estimation of a can 
be drawn from the following considerations. If the 
distance from the adatom to top first neighbours dtop is 
2.6 A < dtop < 3.1 A [231137], the distance to the hollow 
first neighbours is 3.0 A < dhoUow < 

3.4 A. Which 

corresponds to a distance increase of 10 to 15%. The 
overlap integral between s and p orbitals evolving like 
SH, we get 0.75 < a < 0.85. 

In order to gain more insight into how the positions of 
the resonances depend on the model parameters a, ujaa 
and V, figure [5] (c) shows Re[S{uj)] for the top and hol- 
low configurations, for various a > 0.5 and V = 1 eV. By 
graphically solving equation [8] for 1 eV< ujaa < 2 eV, we 



find that increasing a or ujaa increases the resonances 
energy separation and increasing ujaa shifts both reso- 
nances to higher energies. Finally, y is a scale parameter 
appearing as in i?e[S(a;)] (and Im[T,{uj)]) so that an 
increase in V would dilate the vertical scale of the graph 
in figure [5] (c) by V'^. 

We now focus on the residual DOS on the adatom away 
from the resonance at low energy. From equation [7j in 
this energy range, the DOS behavior is determined by 
/m[E(cij)], which is represented in figure [5] (d) for the 
top and hollow configurations, for various a > 0.5 and 
V — 1 eV. Note that Im[T,] is usually interpreted as the 
broadening of the impurity level due to the coupling to 
the continuum. In our case, the situation is more com- 
plex than in the case of a continuum of constant DOS 
with energy as the exotic DOS of graphene introduces 
an energy dependence of Jm[I](w)]. Thus, referring to 
Im[T,(ui)] as the broadening of the adatom resonance is 
an abuse of terminology and we will therefore use quota- 
tion marks. 

Previous studies compared the a — top configura- 
tion case (dotted line) to the a — 1 hollow configuration 
case (black dashed line) and extracted an w dependence 
of /to[I]] in the former case and an oj^/t^ dependence 
in the latter case [16] . The anomalously small "broaden- 
ing" of the adatom band at low energy in the hollow con- 
figuration is explained by the appearance of destructive 
interferences between the different hopping paths con- 
necting three graphene sites of the same sublattice to 
the adatom site (see figure [5] (a)). The adatom in the 
a = top configuration impacting only one graphene 
site, such interferences do not take place and the band 
"broadening" is thus larger [TBHl9j . Introducing hop- 
ping from the adatom to the second graphene neighbors 
in the top configuration (a ^ 0) leads to the same hop- 
ping path symmetry as in the hollow configuration for 
the three second neighbors. Consequently, the same in- 
terference phenomenon happens and contributions up to 
uj^/t'^ appear in the expression of Im[T,'^^°'^^{uj)]. As 
graphically shown by figure [5] (d), the significance of the 
higher order terms depends on the strength of the cou- 
pling to higher order neighbors (aV). Increasing alpha 
decreases (increases) the "broadening" of the top (hol- 
low) resonance. The inset graph shows that the ratio 
Im[I]'^ (lu)]/ Im[I]^ (uj)] diverges at low energy for all con- 
sidered values of a. Hence, we expect that our model 
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FIG. 5: (Color online) The graphene/adatom interaction in the framework of the Anderson model, (a) The top (left) and 
hollow (right) configurations, (b) Interacting level energy in both hollow (uj^j^) and top (ljaa) configurations and their relative 
shift ijJaa ~ ^AA^ for i^aa = 1 eV and V — 1 eV as a function of a. (c) Real part of the adatom self-energy in the top (T^) 
and hollow (Ha) configurations for various a > 0.5 (and a = in the top configuration) and V = 1 eV. The gray straight lines 
represent uj — ljaa for ojaa = 1 eV, 1.6 eV and 2 eV. (d) Imaginary part of the adatom self-energy in the top (Ta) and hollow 
(Ha) configurations for various a > 0.5 (and a = in the top configuration) and V = 1 eV. Insert: Jm[S"^(cj)]/7m[E^(cj)] for 
various a > 0.5 (with the same color code as for /m[S]). (e) Normalized DOS associated to the adatom in both top and hollow 
configurations with ljaa ~ 1.6 eV, V = 1 eV and a = 0.8. Insert: Residual DOS at low energy. 



accounts for the larger residual DOS for the top adatom 
than for the hollow adatom at sufficiently low energy. 

Figure [5] (e) finally shows the density of states as- 
sociated to the adatom for both configurations, with 
^AA — 1-6 eV and a — 0.8. These values have been 
obtained by solving equation |8] for the hollow and top 
configurations with = 1.35 eV and = 1.6 eV 
in order to fit the experiment (see figure Hid)). The in- 
sert displays the calculated DOS at low energy. In this 
energy range, the DOS is very low with respect to the 



vicinity of the resonance, which is consistent with the 
typical evolution of the tip-surface distance as a function 
of tunneling bias reported in figure S2. When consid- 
ering the DOS ratio, we find that the DOS on the top 
adatom clearly dominates the one on the hollow adatom 
(at UJ = 0.6 eV, p^^/pS « 7) [42]. Thus, three different 
energy domains emerge, successively where the DOS of 
the top adatom dominates (oj < 1.0 eV in this case), then 
the DOS from the hollow adatom gains the upper hand 
(1.0 eV< UJ < 1.5 eV) and again the DOS from the top 
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adatom dominates (w > 1.5 eV). This result is thus con- 
sistent with the experimental measurements presented in 
figure |3] showing two moire contrast inversions. In the 
end, with reasonable values for a, V and ujaa, the simple 
model we used accounts for the variations in level energy 
and broadening for the two different graphene/adatom 
stacking configurations expected from the experiments. 

To finish, note that by comparing the calculated LDOS 
to the experimental conductance spectra, we neglect the 
energy-dependence of the tunnel barrier transmission co- 
efficient T{E, Vs, z) whose effect is to favor tunneling to 
higher energy empty states i.e. close to the resonances. 
Also, we assume a weak coupling to the tip states. More- 
over, we do not take into account the Fano-resonance 
effect that can arise from interference between tunnel- 
ing paths connecting the tip and graphene and the tip 
and the localized state [18j [191 ES]- In fact, we expect 
no such effect here because graphene on SiC is electron- 
ically transparent on high bias STM images ^28, t44» i45j . 
Indeed, the STM constant current images of figure |3] do 
not show the graphene atomic contrast for tunneling bi- 
ases down to Vs — -1-0.5 eV while the adatom states are 
clearly identified. In other words, at high tunneling bias, 
direct tunneling from the tip states to the adatom states 
is highly favored with respect to tunneling from the tip 
states to the graphene states even though the adatoms 
lay below the graphene plane. 



D. The graphene-adatom interaction from ab initio 
calculations for graphene on SiC(0001)(2 x 2)c 

In this section, we study the graphene-adatom inter- 
action using ab initio calculations that take into account 
the complete graphene-SiC(0001)(2 x 2)c system. A de- 
tailed description of the calculation was published in [21] . 
Here we focus on the energy range that corresponds to 
the graphene-adatom interaction. Due to technical con- 
straints on the calculation supercell size, we chose the 
quasi-commensurate 5x5 graphene on top of a 4 x 4 
SiC cell. In order to examine various graphene/adatom 
stacking, we consider two supercells which correspond to 
translations of the graphene lattice with respect to the 
substrate, as shown in figure |6] (a) and (b). With re- 
spect to the graphene lattice, the A4 adatom sits in a 
top configuration, the Al-3 adatoms in a sligthly shifted 
top configuration, the B4 adatom in a hollow configura- 
tion and the Bl-3 adatoms in a bridge configuration. 

The supercells we consider correspond to two specific 
geometries, each one giving rise to a small moire pat- 
tern. It is different from the geometry probed by STM 
which gives rise to a large moire pattern that contains all 
the graphene-adatom stacking configurations listed above 
within one pseudocell. We thus use our calculations to 
study the graphene-adatom interaction rather than as a 
model of the system probed experimentally. Moreover, 
we have to keep in mind that our GGA calculations of 
the graphene-SiC(0001)(2 x 2)c may have some imper- 



fections arising from the possible underestimation of the 
graphene-reconstruction interaction [37], which however 
are not critical for the purpose of this study. Conse- 
quences of this are the following. Firstly, we cannot dis- 
cuss the graphene-restatom interaction since our calcu- 
lations do not give account for this small coupling [24j . 
The LDA and GGA(D) calculations from the literature 
show a band anticrossing of ~ 50 meV [23 and we 
have shown experimentally that the graphene-restatom 
interaction indeed exists in the discussion in connection 
with figure [3] However, we chose to concentrate on the 
much stronger graphene-adatom interaction, which we 
can discuss with our GGA calculations. Secondly, we ob- 
tain a neutral graphene layer [24] while angle resolved 
photoemission spectroscopy [55* and calculations cited 
above 37 indicate electron doping of the graphene sheet 
so that Ep ^ Ed+ 0.2 eV and Ep ^ Ed + 0.35 eV re- 
spectively. Thus, our calculations may not reproduce the 
position of Ep in this system. Finally, we obtain a flatter 
graphene layer (0.07 A [23]) than calculations [37] from 
the literature which show a small corrugation of 0.16 A 
using LDA and 0.27 A using GGA(D). We will discuss 
the consequences of this later in this section. 

The calculated band structure for each configuration 
is given in figure [6] (c) and (d). In both cases, it exhibits 
a clear band anticrossing between the 4 weakly disper- 
sive adatom bands and the graphene tt* band, evidencing 
interaction between graphene and adatom states. Note 
that far from the adatom states, near the Dirac point, 
our calculations in the same system showed that the lin- 
earity of the TT bands is recovered [53] , which is also the 
case for GGA(D) and LDA calculations '37'. 

Figure[6](e) gives the calculated adatom LDOS for each 
type of stacking. Features displayed by figure |6] (e) are 
complex but the simple model introduced in the previous 
section will help analyze them. By comparing the two ex- 
treme situations, top stacking (red) and hollow stacking 
(black), we find a lower position in energy of the spec- 
tral weight for the hollow than for the top configuration. 
Regarding the broadening of the adatom band, we find 
a larger extension in energy in case of the top configu- 
ration. Hence, we find the same trends in the impact of 
the graphene/adatom interaction on the adatom level as 
with the Anderson single-impurity model. 

Considering intermediate configurations which we did 
not treat with the previous model, the slightly shifted 
top position leads to a LDOS very similar to the one for 
the exact top configuration, in energy and band broad- 
ening. Such behavior proves that the graphene/adatom 
interaction is not strictly local. The bridge configura- 
tion leads to an average band energy position similar to 
the the one for the top configuration but with a slightly 
narrower distribution. 

Going more into details, the fine structure of the ab 
initio calculated LDOS may arise from various factors 
that were not taken into account in the previous model. 
Indeed, in this system we consider an adatom array on 
a semi-infinite medium, while in the previous section the 
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FIG. 6: (Color online) Ab initio electronic structure cal- 
culations of graphene on SiC(0001)(2 x 2). (a) and (b) 
Top view of the two considered supercells, displaying top 
(A4), slightly shifted top (Al-3), bridge (Bl-3) and hol- 
low (B4) graphene/adatom stacking configurations. Only 
the last SiC bilayer, the adatoms and the graphene layer 
are represented, (c) and (d) respective band structure cal- 
culations in the energy range spanning the adatom low- 
dispersive bands, (e) Calculated LDOS for the A3-4 and 
B3-4 adatoms. (f), (g) Maps of l^^l^ integrated over the en- 
ergy range [-1-0.5 eV, -1-1.0 eV] taken just above the graphene 
atoms, for the configurations given in (a) and (b). Contrast 
is reversed with respect to Fig. |2] (d) and|3](b): higher den- 
sity is darker. The adatom positions are indicated by dashed 
circles. 



impurity was isolated. The weak direct adatom/adatom 
interaction can slightly affect the results. Moreover, each 
adatom interacts with graphene which is in interaction 
with the other adatoms. In the end this can result in 
the multiplication of the peaks in the adatom LDOS by 
supercell (periodicity) effects. Since the multiplication of 
the peaks is governed by the calculation geometry, which 
is different from the one in the experimentally probed 
system, we cannot reasonably extract any general infor- 
mation from the fine structure of the calculated LDOS 
spectra. Moreover, STM measurements were conducted 



at room temperature, which implies a 100 meV broad- 
ening in energy of the LDOS features. For instance, the 
LDOS crossing at 0.69 eV in [6] (e) is unlikely to appear 
as a contrast inversion in our STM measurements. 

Finally, figures [g] (f-g) show maps of j^E'p integrated 
over the energy range [-1-0.5 eV, -1-1.0 eV], taken just 
above the graphene atoms, for the configurations given 
in (a) and (b). We chose the energy range where figures 
|6] (c-d) show a clear hybridization between the graphene 
and adatom states, in order to reveal the impact of the 
graphene-adatom interaction on the graphene electronic 
structure in the direct space and compare it to STM data. 
Note however that the integration energy range in the ab 
initio calculated maps is different than the one on low 
bias STM images. Firstly, it is difficult to compare the 
very low integrated DOS calculated for lower energy with 
respect to the Dirac point to STM data. Secondly, we re- 
call that in the experimental system, the position of Ep 
is most likely shifted toward higher energy with respect 
to the GGA calculation. 

Figures [6] (f-g) show LDOS reduction on the graphene 
atom atop of the adatom (figure [6] (f)) and equivalent 
LDOS on the 6 neighboring graphene atoms in the hol- 
low configuration (figure |6] (g)). We thus find the same 
kind of features we identified on low bias STM images 
(figure [2] (d) and[3](b)). Going more into details, on the 
low bias STM images, the dark contrast not only corre- 
sponds to top, but also to side top and to bridge sites 
of the Si adatoms. The shape and apparent depth of 
the "switched off" atoms depends on the local stacking. 
Only for nearly exact hollow site configuration does the 
perturbation disappear. First, this is consistent with the 
fact that the graphene-adatom interaction has a large 
spatial extension and thus we consider a ~ 1 in the An- 
derson model. It is also consistent with the fact that 
adatoms in top, side top and bridge configurations are 
similarly impacted, as indicated by our ab initio study 
of the adatom LDOS (figure |6] (e)). Moreover, we recall 
that in the Anderson model, it is only when the very spe- 
cific hollow configuration is reached that the interference 
between the hopping paths (at low energy) arises, and 
that the graphene-adatom states hybridization is subse- 
quently strongly reduced. 

To finish, we pointed out that our calculations possi- 
bly underestimate the graphene-reconstruction interac- 
tion and consequently underestimates the graphene cor- 
rugation. Our results thus show that the extinction of 
the LDOS on the graphene atom in top configuration is 
already present for a fiat graphene layer. For a corru- 
gated surface, the top graphene atom is closer to the Si 
adatom by less than 0.5 A according to our low bias mea- 
surements and ab initio calculations [37| . This translates 
in a larger value of V in the Anderson model i.e. an 
even larger graphene-adatom states hybridization. From 
the distance dependence of tight binding parameters, V 
should however increase by less than 25%. We recall that 
STM constant current images represent both the topog- 
raphy and the electronic structure of the surface. In the 
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end, we expect an enhanced contrast on low-bias STM 
images from the topographic lowering of the top graphene 
atoms. Regarding the Si adatom resonances, a 25% larger 
V in the top configuration reduces the hollow-top reso- 
nance energy difference - but not enough to invert their 
positions (see figure [s] (c))- and further increases the top 
resonance "broadening" with respect to the hollow res- 
onance (see figure [5](d)). Thus, taking the corrugation 
into account will not strongly modify the conclusions of 
this study. 

IV. CONCLUSION 

We have been able to investigate the graphene- 
impurity interaction as a function of the graphene- 
impurity local stacking. A direct experimental com- 
parative study has been possible thanks to the speci- 
ficities of graphene on SiC(000T)(2 x 2)c. The sub- 
strate reconstruction is composed of one Si adatom 
and one C restatom per unit cell, and is mainly in- 
teracting with graphene through the adatoms. It has 
a large cell parameter of 6.14 A, allowing us to con- 
sider the adatoms as independent in a first approxi- 
mation. The graphene/adatom interaction is apprecia- 
ble but not too strong so that the general electronic 
structure of graphene and the reconstruction are main- 
tained. The rotation between the graphene and the re- 
construction lattices imply spatial variations of the local 
graphene/adatom stacking, allowing us to compare dif- 
ferent configurations on the same STM image. Finally, 
the electronic transparency of graphene on SiC on high 
bias STM images makes it possible to clearly identify the 
local stacking and to study both the impact of the inter- 
action on the graphene LDOS - low bias tunneling regime 
- and on the adatom level - high bias tunneling regime. 

Low bias STM images show that the graphene LDOS 
near the Fermi level is clearly affected in the top- 
like graphene/adatom stacking regions, graphene atoms 
atop of an adatom showing strongly reduced LDOS, 
while perturbations are hardly noticeable in hollow 
graphene/adatom stacking regions. Regarding the 



adatom resonance, a combination of constant current 
STM images at various high tunneling biases and tunnel- 
ing spectroscopy measurements reveals a lower resonance 
energy and smaller broadening in the hollow configura- 
tion than in the top configuration. 

The single-impurity Anderson model gives consis- 
tent results with the experimental observations for the 
adatom resonances in the top and hollow configurations 
when the same interaction cut off for both configurations 
is considered. Within this model, the reduced impu- 
rity level broadening in the hollow configuration arises 
from interferences between hopping paths connecting the 
adatom site and the neighboring graphene atoms within 
each graphene sublattice. 

Finally, we theoretically investigate the graphene- 
adatom interaction using ab initio calculations that take 
the whole graphene/reconstruction structure into ac- 
count. Considering the LDOS on adatoms, from the 
more complex calculated electronic structure, we extract 
the same trends in the stacking-dependent impact of the 
interaction on the adatom level as with the simple An- 
derson model and STM. 

To conclude, our results confirm the importance of lo- 
cal stacking in the graphene-impurity problem. They 
show that it is possible to discriminate between two lo- 
cal stacking configurations by comparing the impurity 
level broadening and shifting. The local stacking can 
also be identified by considering the perturbations in the 
graphene LDOS. 
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